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Abstract 

Higher moments of the vorticity field ft m (t) in the form of L 2m -norms (1 < m < oo) 
are used to explore the regularity problem for solutions of the three-dimensional incompressible 
Navier-Stokes equations on the domain [0, L]^ er . It is found that the set of quantities 

1m 

provide a natural scaling in the problem resulting in a bounded set of time averages (D m ) T on a 
finite interval of time [0, T]. The behaviour of D m +i/ D m is studied on what are called 'good' 
and 'bad' intervals of [0, T] which are interspersed with junction points (neutral) tj. For large 
but finite values of m with large initial data (O m (0) < woO(Gr 4 )) , it is found that there is an 
upper bound 

< cl v zu a Gr 4 , zu = i>L~ 2 , 

which is punctured by infinitesimal gaps or windows in the vertical walls between the good/bad 
intervals through which solutions may escape. While this result is consistent with that of Leray 
[lj and Scheffer [10J, this estimate for O m corresponds to a length scale well below the validity 
of the Navier-Stokes equations. 
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1 Introduction 

The challenge that analysts have faced in the last 75 years has been to prove the existence and 
uniqueness of the three-dimensional Navier-Stokes equations for arbitrarily long times [1] [21 OH EE]. 
Its inclusion in the AMS Millenium Clay Prize list [7j has widely advertised the the nature of the 
problem but the elusiveness of a rigorous proo{j] and the severe resolution difficulties encountered in 
CFD, even at modest Reynolds numbers, are puzzles that have grown as the years progress. 

Nevertheless, there is a long-standing belief in many scientific quarters, on the level of a folk- 
theorem, that the three-dimensional Navier-Stokes equations 'must' be regular. Mathematicians 
are more cautious and still take seriously the possibility that singularities may occur, at least in 
principle. Leray p] and Scheffer [10] proved that the (potentially) singular set in time has zero 
half-dimensional Hausdorff measure [IT]. The Leray-Scheffer result motivated Caffarelli, Kohn and 
Nirenberg [12J to introduce the idea of suitable weak solutions to study the singular set in space-time 
which they concluded has zero one-dimensional Hausdorff measure. Thus, if space-time singularities 
exist then they must be relatively rare events. These ideas have spawned a growing literature on 
the subject where more efficient routes to the construction of suitable weak solutions are in evidence 

It is worth remarking that the wider issue regarding the formation of singularities has been 
obscured by the very great difficulty that exists in distinguishing them from rough intermittent 
data. Intermittency is characterized by violent surges or bursts away from averages in the energy 
dissipation, resulting in the spiky data that is now recognized as a classic hallmark of turbulence 
[22j[23l[2l[25j[26]. At least three options are possible: 

a) Solutions are always smooth with only mild excursions away from space and time averages; 

b) Solutions are intermittent but, despite their apparent spikiness, remain smooth for arbitrarily 
long times when examined at very small scales; 

c) Solutions are intermittent but spikes may be the manifestation of true singularities. 

Options b) and c) are impossible to distinguish using known computational methods. The Leray- 
Scheffer result shows that potential singularities in time must be distributed as no more than points on 
the time axis, but it contains little other information. Both for analytical and computational reasons 
it would be desirable to understand the origin of these points and the structure of the solution near 
to them. The aim of this paper is to address this issue. 

In the past generation physicists have used Kolmogorov's theory to examine intermittent events 
by studying anomalies in the scaling of velocity structure functions. This theory is based on a set of 
statistical axioms, not directly on the Navier-Stokes equations. Nevertheless, to make a comparison, 
the intermittent dynamics discussed above would lie deep in the dissipation range of the energy 
spectrum. Frisch's book [27] and the recent review by Boffetta, Mazzino and Vulpiani [28] contain 
readable accounts of these ideas. 

1 Cao and Titi [8] and Kobelkov [9] have recently proved the regularity of the primitive equations of the atmosphere 
and oceans, even though these have been considered by many to be a problem harder than the Navier-Stokes equations. 
The methods used unfortunately do not appear to successfully transfer to the Navier-Stokes equations. 
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1.1 General strategy 

The main idea of this paper is to use higher moments of the vorticity field u> instead of derivatives. 
Scaled by a system volume term L -3 , a set of moments with the dimension of a frequency are defined 
such that for m > 1 

n m (t) = lL- 3 j^\u\ 2m dV\ +w , (1.1) 

where zuq = uL~ 2 is the basic frequency of the domain of side L. Q,\ is synonymous with the 
i^i-norm and sits within the sequence of inequalities 

w Q < n x (t) < n 2 (t) <■■■< n m (t) < n m+1 (t) <..., (1.2) 

so control from above over Q m for any value of m > 1 implies control over the i/i-norm which, in 
turn, controls from above all derivatives of the velocity field [2| [3l [5l [6] . 

A technical problem lies in how to differentiate the ^l m (t) and manipulate them without the 
existence of strong solutions for arbitrarily large t. This difficulty can be circumvented my restricting 
estimates to a finite interval of time [0,T] and then pursuing a contradiction proof in the following 
standard manner. Assume that there exists a maximal interval of time [0,T max ) on which solutions 
exist and are unique; that is, strong solutions are assumed to exist in this interval. If [0,T max ) is 
indeed maximal then Qi(T max ) = oo. The ultimate aim of such a calculation would then be to show 
that lim sup T ^ Tmaa: ^ m ' s f' mte f° r anv m> 1; \f this turned out to be the case it would lead to a 
contradiction because [0, T max ) would not be maximal. Thus T max must either be zero or infinity: it 
cannot be zero because it is known that there exists a short interval [0, irj) on which strong solutions 
exist, so T max = oo. 

The results in Sj2]have been estimated using this strategy. It turns out that there exists a natural 
scaling within the Navier-Stokes equations which makes the variable 

D m (t) = (wQ l 9. m ) am with a m = 2m , (1.3) 

the most natural to choose. Then Theorem Q] shows that 

(D m ) T <c av Gr 2 + 0{T- 1 ), (1.4) 

with a uniform constant c av . Two remarks are in order. Firstly it is not difficult to extract an estimate 
for a set of length scales from (Oil . Defining A~ 2c * m = v~ a,n (f^ m ) T , this shows that 

{ L Kn ) 20tm = (Dm) T , (1-5) 

and thereforaEI 

LX- 1 < (c av Gr 2 ) 1/2am + OiT- 1 ) . (1.6) 

The exponent a m within the definition of D m appears to be a natural scaling consistent with that 
of the Sobolev inequalities. This paper suggests that the breaking of this scaling through stretching 



2 Doering and Foias [37j have shown that for Navier-Stokes solutions Gr < cRe 2 which would be valid if solutions 
were assumed to exist for large enough values of T. In this case the Gr 2 -term on the right hand side of (|1.4p would 
be replaced by Re 3 in which case the right hand side of (| 1 . 6j) would be Re 3 ^ 2am . Thus, LAj" 1 < c 1 ^ 4 Re 3 ' 14 ' which is 
the Kolmogorov estimate. For large m, this becomes significantly larger running to LA^ 1 < cRe 3 . 
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between D m +i and D m may be required to make progress. This is gauged more specifically in 
Theorem [2] in where it is shown that a finite interval [0, T] of the time axis can be potentially 
broken down into three classes, denoted by good and bad intervals with set of junction points (or 
intervals) {r^} designated as neutral. In it is found that the direction of the inequality is reversed 
on the good and bad intervals; that is 



Dm+1 < c D -»mQ r P{T) 

Urn 



< (good) 

= (neutral) (1.7) 
> (bad) 



In (II. 7ft p(T) is a T-dependent exponent (> 2) of the Grashof number Gr and p m is a parameter 
in the range < fj, m < 1. The universal inequality £l m < fi m +i ultimately shows that on good and 
neutral intervals 

D m < Q m m , (1-8) 

where Q m is a function of p(T), Gr, a m and fi m . The main question lies in the nature of the 
transition from the good to the bad intervals through the neutral points 73. On bad intervals the 
application of the reverse inequality in (jl. 71) to the differential inequality for D m in Proposition Q] 
results in regions smaller in amplitude than Q m in which solution trajectories remain bounded by 

D m < B m ™ . (1.9) 

The bad regions are not absorbing : solutions remain inside these regions if they enter inside, but they 
are not attracted into them if they lie outside. The key point is that for all finite values of m > 1, 
£>m < Gm, thereby leaving vertical gaps or windows through which trajectories can potentially escape 
to infinity - see Figures 1, 2 and 3. However, while the gap between Q m and B m closes for large 
m, the limit m = 00 is forbidden and so these windows can only be reduced to infinitesimally 
small holes which puncture a general upper bound. This result is consistent with that of Leray 
[lj and Scheffer [10]. In terms of Cl m , this punctured bound turns out to be 

On < c 2 av w Gr 4 (1.10) 

When converted into a length scale, this estimate shows that regular solutions may go as deep as 
near nuclear scales (10~ 2 angstroms) and therefore many orders of magnitude below the validity of 
the Navier-Stokes equations. The conclusion is that unless other unknown controlling mechanisms 
are shown to exist, the Navier-Stokes equations may formally possess solutions that either become 
singular or, if they continue to exist, may be unresolvable numerically. 

1.2 Notation and functional setting 

The setting is the incompressible (divit = 0), forced, three-dimensional Navier-Stokes equations for 
the velocity field u(x, t) 

u t + u ■ V« = uAu - Vp + f(x) , (1.11) 
with the equation for the vorticity expressed as 



u: t + u ■ Vu; = vAu + u ■ Via + curl/ . 



(1.12) 
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The properties of the forcing & other definitions are given in Table [TJ The domain V = [0, L] 3 
is taken to be three dimensional and periodic. The forcing function f(x) is L 2 -bounded and the 
Grashof number Gr is proportional to ||/||2 : see the paper by Doering and Foias [37] for a discussion 
of narrow-band forcing [37J : for simplicity the forcing is taken at a single length-scale £ = L/2tt. 



Quantity 


Definition 


Remarks 


Box length 


L 




Forcing length scale 


e 


t = L/2vr 


Average forcing 


frms = L 3 11/ III 




Narrow-band forcing 


ii/iii«mv n /ui 


n > 1 


Grashof No 


Gr = ef rms v- 2 




Box frequency 


wq = vL~ 2 




Characteristic velocity 






^-definition 


E(t) = f v \u\ 2 dV 


Energy 


/3 m -definition 


(3 m = m(m + 1) 




a m -definition 


— 2m 
«m — 4m _ 3 




p m -definition 


p m = 2m(4m + l)/3 





Table 1: Definitions of the main parameters. The forcing is taken at a single length-scale . 
Now define 



\uj\ 2m dV. 



where the frequencies fi m are given by 

n m (t) = (l~ 3 j m ) 1/2m + m . 

The term wq in (j!.13p provides a lower bound for Q m . Indeed it is easy to prove that 

^0 < ni(t) < n 2 (t) <...< n m (t) < n m+1 {t) <... 

The symbol ( • ) T denotes the time average up to time T 

cT 



1 r 

(9(-))t = limsup- / g{r)dr. 
9(0) 1 Jo 



L/2n. 
(1.13) 
(1.14) 
(1.15) 

(1.16) 



9(0) 

2 Some properties of the fl m (t) 

2.1 A differential inequality and a time average 

This subsection firstly contains a result concerning the differential inequalities that govern the set of 
frequencies Cl m (t). Secondly it contains a result that is an estimate for an upper bound on a set of 
time averages over the interval [0, T]. Finally it contains a result on the nature of exponential bounds 
on [0, T]. All of the proofs, which lie in Appendices A, B and C, are based on the contradiction 
strategy explained in §1.11 Firstly we define 

2m 

a m = - -. (2.1) 



Dn 



(w fi r 



4m — 3 
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Proposition 1 On [0,T], for 1 < m < oo, n = \{m + 1) and Gr > 1, the D m satisfy 



(w a m ) D m < D Ti 



1 / D 



J m+1 



Dl, + c 2 , m Dl + c 3 , m Gr } , 



(2.2) 



where p m = 2m(4m + l)/3. For the unforced case the last term on the right hand side of \2. 2|) is 
proportional to c^ m . 

Remark: Note the strict inequality m < oo : the Riesz transform used in the proof in Appendix [Al 
requires the introduction of higher derivatives when m = oo. 



Theorem 1 : For 1 < m < oo and Gr > 1 



(D m ) T < c av ( Gr 2 H 3 — - 



where Eq = E(0) is the initial value of the energy. For the unforced case, the estimate is 



{D m } T < C 



L~ 5 E 
wlT 



(2.3) 



(2.4) 



Remark: (12. 3p can also be expressed as 

(D m ) T < c av Gr p , (2.5) 
where C is a uniform constant. The m-independent exponent p(T, Eo,Gr) written as 

p(T, E , Gr) = 2 + In (l + ^J^ Gr ~A (inGr)- 1 . (2.6) 

3 Trajectories on good, bad and neutral intervals 
3.1 The ratio D rn+1 /D m 

Given the result in Proposition [TJ understanding the behaviour of the ratio D m j 1 D m+ \ is an important 
step. 

Theorem 2 For the parameters p m = p, m (T, p, Gr) with values in the range < p m < 1, the ratio 
D m /D m+ x obeys the inequality 



Dr. 



D 



m+l 



l Gr- p( ^D!i" 



(1— flm)//Jm 



> 0. 



(3.1) 



Remark 1 : The proof lies in Appendix[C] and is dependent on the result of Theorem [TJ 

Remark 2: Theorem [2] implies that while there must be intervals where the integrand is positive, 
there could also be intervals where it is negative. While it tells us nothing about the interval size or 
distribution it is clear that these are T-dependent. 
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Formally the theorem leads to the conclusion that there exists at least one good interval of 
time within [0, T] on which : 



D r 



D 



> 



m+l 



-i 



while there potentially exist bad intervals of time on which 



D r 



D r 



- < \c av Gr^ 
•l L 



-i 



Neutral points or intervals represented^ by the zeros of the integrand in (13. ip lying at 

n = Ti{n m , P(T), Gr) . 
In terms of ^l m +i and Q m (13.21) and p. 31) become 







a, 



where Q m and 7 m are defined by 



m+l < 



> 



good 

neutral 

bad 



3 

Pm — l^m 



m(4m + 1) 

The positivity of 7™ requires that /x m be bounded away from zero such that 

3 

< H m < 1 . 



m(4m + 1) 

Beca use f2 m +i > ^m> (j3.5|) shows that on good and neutral intervals 

D m ,good<GT l<m<(XD. 

Now we turn to the bad intervals : consider (13. 3p in (12. 2p . in which case (Sl n < U m ) 



{zu a r 



(3.2) 
(3.3) 
(3.4) 

(3.5) 

(3.6) 
(3.7) 
(3.8) 

(3.9) 

(3.10) 



{-^- (c av Gr^D7 n f™) Pm D 2 m + c^D 2 ^*™ + c 3 , m Gr| , (3.11) 



where p m = 2m{Am + l)/3 but m = oo is forbidden. The range of validity of [i m in (j3.9j) can be 
re-written as p m > /i m/ o m > 2. Thus -D m < if, at the time of entry n into a bad interval 



c av Gr^D mMd { n ) 



D m > Ci im C2,mD m ,bad(Ti) 



(3.12) 



Given that p m [i m > 2 and a n > a m , the first term on the right hand side of (j3. 12H is dominant. 
Using the lower bound D m > 1 it is found that 



D 



m, bad 



where 



B n 



1 



^Gr^ 

2(dn Gro) ~\~ OirnPm^i 



-Gr 



l/a„ 



Cl,m,C2,m 

bm — O-mPmPm 



(3.13) 

(3.14) 
(3.15) 



There is no information on how the t% are distributed. 
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3.2 How large are Q\ and B\l 

Di(t) 
Gr 5 



Gr 2 



i 

1 




good 

s / 




\ good 




\ / 




\ 




•/•V- j 











bad 




bad 



t 

n r 2 r 3 

Figure 1 : From a variety of initial conditions for m = 1 t/?e cartoon above shows how solutions may 
potentially escape at or near neutral points t = t\ or a later value t = T3, or even return att = t%. However, 
all must satisfy the bound on the time-average. 

For m = 1 we have b\/a\ = 1 and p\ = 10/3 ; the difference in the sizes of Q\ and B\ lies in the 
upper bounds on [i\ and on The latter has been defined in (|3.8p 



Hi < 1 



Hi < 1 -3/5 = 2/5. 



From (i3~BT) and (}3~l0j) we have 



Di, good < (cavGrV) 1 '^ 
which, on minimization of the right hand side, gives 

D l>good < {c av Gr 2 f 2 = cH 2 Gr 5 . 

The equivalent estimate for Di^ad is 



D 



1, fiad 



< 



( Cljl c 2jl ) 3 A° 



Gr 2 - 0(Gr 



3/10n 



(3.16) 
(3.17) 

(3.18) 
(3.19) 



It is useful to re-work these estimates in terms of a point-wise invers43 length-scale ■q l 4 = v 3 e with 
a point-wise energy dissipation rate e = uQf = vw\D\. The result, 

Lvi 1 < c^Gr^ (3.20) 
is shown in Figure 3 where the constant on the bad estimate is slightly smaller. 



4 The context of this is the estimate for the inverse length LA^ < cJ^Gr 1 ^ of JT] 
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Lvl 1 








Gr 5 / 4 








good 


good 




Gr 1 ' 2 




(time average) 








bad 


bad 



n t 2 r 3 

Figure 2 : Bounds on Li]^ 1 : notice the large size of the gaps between the good and bad intervals. Based 
on the constants, the upper bound on the time average is larger than that on the bad intervals. 



3.3 How large are and B% n for large ml 



D m (t) 



Gr 2 



small gaps through which 
trajectories may pass 



1 t 

Tl T 2 T 3 

Figure 3 : For large m, the gap between C?," m and is infinitesimally small but the limit m = oo is forbid- 
den. The upper bound on the time-average is the horizontal line of dots. At t\ and T3 a solution must enter 
the corresponding bad interval within the upper bound to remain inside. 

From the definitions of (j3.6j) and ([3.140 and the fact that jl m < fj, m , it is clear that Q^j 1 —6%™ > 0, 
keeping in mind that the limit m = 00 is forbidden. Firstly the Cj jTO are polynomial in m and 
Pm ~ 0(m 2 ) for large m. Therefore 

(ci,mC2,m) /* 1 , b m /a m /l, and fj, m / Ji m . (3.21) 
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Hence, for large m 

g%>-B%>\0. (3.22) 

Specifically for D m , for very large m, the upper bounds on fi m and fl m can now be taken arbitrarily 
close to unity provided that fi m < 1 and p, m < 1 . From (j3.6|) , minimization of the right hand side 
gives 

D m , good <c av Gr 2 (\). (3.23) 

The equivalent estimate for Di^ad is 

An, bad < ° av , Gr 2 / c av Gr 2 . (3.24) 

(Cm,lCm,2) /Pm 

4 Conclusion : what are the length scales corresponding to the upper bounds? 

The key feature of this paper is the closure of the gaps between the good/bad intervals as m —> oo 
but with the actual limit m = oo forbidden. The origin of this lies in Proposition [TJ in the use of the 
inequality (p = \{m + 1)) 

||Vu|| p < Cp||u;||p pe(l, oo), (4.1) 

whereas, when m = oo 

UVulloo < cllwIlooCl + lnfla) • (4.2) 

(14. ID has its origin in a double Riesz transform while (14. 21) arises from the work of Beale, Kato and 
Majda [38] on the three-dimensional Euler equations - see also Kato and Ponce [39]. The InH^ term 
in (j4.2j) prevents the closure of the set of inequalities for D m . While the m = oo limit is valid for 
good intervals, it is not valid for the bad because of the necessary use of Proposition [TJ Thus it is not 
possible to completely close the gaps between the two sets of intervals, although they can become 
arbitrarily small. This allows for the possibility of the escape of trajectories. The m-dependence of 
the Ti means that the junction points can, in principle, lie at different places on the time-axis as m 
varies. If the gaps fall randomly with respect to m then a trajectory would have to thread its way 
through these to escape to infinity. However, an unknown but subtle alignment of the gaps cannot 
entirely be ruled out. 

The closeness of the upper-bounds on both the time average and on point-wise values of D m 
(m S> 1) away from the gaps, poses the question whether there exists dynamics that naturally lie 
either close to these bounds or even fulfill them. The point-wise energy dissipation rate per unit 
volume is 

£ = vVi 2 < v 3 L- 4 dV<*™ u ^L- A c av Gr 8 . (4.3) 
Defining a local Kolmogorov length as \kjoc = {z/v 3 ) 1 ^ we obtain 

L\^ loc < c av Gr 2 , (4.4) 

which is consistent with the estimate in (jl.5p for large m. If the solution survives for large enough T 
to make sense of a Reynolds number based on Uq = L~ 3 (\\u\\?,) T , then the Doering-Foias result for 
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Navier-Stokes solutions [37], Gr < cRe 2 , can be invoked to give an estimate for a local Kolmorgorov 
scale@ 

L\ k j oc <cRe 3 . (4.5) 

In the atmosphere, for instance, this length-scale would be of O(10 -12 ) metres - about 10~ 2 angtroms 
- which is about the scale of the nucleus (!) and is thus outside the validity of the Navier-Stokes 
equations. Because the bounds on the good and bad intervals are very close to the time average then 
solutions could, in principle spend long periods of time close to this bound and remain regular, yet 
such a scale is not only unreachable computationally but is outside the validity of the NS equations. 
Thus, a singularity is not necessary to produce unresolvable solutions. 

Acknowledgements: I would like to express very warm thanks to Claude Bardos, Matania Benartzi, 
Toti Daskalopoulos, Darryl Holm, Roger Lewandowski, Gustavo Ponce, James Robinson and Edriss 
Titi for discussions on this topic. 



A Proof of Proposition Q] 

Consider the time derivative of J m defined in (II. 13ft 

— j m = ! M 2(m-1) w-{i/Aw + w-Vu + curl/} dV. (A.l) 
2m J v 

Bounds on each of the three constituent parts of (lA.lj) are dealt with in turn, culminating in a 
differential inequality for J m . In what follows, c m is a generic m-dependent constant. 

1) The Laplacian term : Let <p = to 2 = u> ■ u> . Then 

f |w| 2 f m_1) w. AudV = f ^-^w-Aw. 
Jv Jv 



>dV 



f <j) m - l {&{i<j))-\Vu>\ 2 }dV 
Jv 



< I m - x A(i0)dy. (A.2) 



2 

V 



Using the fact that A(</> m ) = m{(m - l)0 m - 2 |V^| 2 + (p^A^} we obtain 

f \u\ 2 ^ m -^u ■ AudV < -\{m-l) [ (j) m - 2 \V<i)\ 2 dV + — f 
Jv Jv 2m Jv 

Jv 



AU m ) dV 



m 2 



2(m - 1) 



m 2 



V{u m )\ z dV. (A.3) 



Thus we have 



5 The correspondence is that Gr 2 is replaced by Re 3 . 
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where 



m 



A m = UJ m C!,m = 7 , (A.5) 

m — 1 

where there is equality for m = 1. The negativity of the right hand side of (1A.4I) is important. Both 
||VA m ||2 and HAnlb will appear later in the proof. 

2) The nonlinear term in (lA.lj) : The second term in (jA.l|) is 

J \u;\ 2( - m - l) u-(u-V)udV < c m (J \u\ 2{m+1) dV^j m+1 (J \Vu\ m+1 dV^J m+1 

m 1 

< c m (j \u\ 2{m+1 UvY +1 (j \u\ m+1 dvY l+1 (A.6) 



where the inequality ||Vu|| p < c p ||cj||p for p G (1, oo) has been useqj: this can be proved in the 
following way: write u = curl(— A) -1 ^. Therefore u» j = RjRitUi where Ri is a Riesz transform. 

Together with (IA.2P this makes (jA.lj) into 



Jm ^ 

2m ci 



/ m 1 r 

\V{u m )\ 2 dV + CmJm+{ J|£ +1) + J \u>\ 2 ^u> • curl/ . (A.7) 



,m JV 



3) The forcing term in ( lA.lj) : Now we use the narrow-band property of the forcing (see the Table 
in £ 11.21) to estimate the last term in (IA.7P 



[ \u\ 2{m -V u -cur\ fdV = [ H 2 ^" 1 ^ -curlew 
Jv Jv 



< 



U 2m dV 



(2ra-l)/2m 



\ l/2m 

|v/| 2m dy] 



(A.8) 



However, by going up to at least 3-derivatives in a Sobolev inequality it can easily be shown that 
||V/||2m < c||/||2-k 2m , because of the narrow-band property. (IA.8jl becomes 



2(m-l) 



a; • curl/ <iy 



2m-l 

< c (L 3 ^) 2 - ||/|| 2 i/ 



— Li 'm Jrms ±J 

< cVL 2 ™- x L z wlGr 



(A.9) 



4) A differential inequality for J m ■ Recalling that A m = to m 

J m+1 = J \u\ 2 ^ dV = J \A m \ 2{m+l)/m dV = \\A m \\X + + %t ■ (A-10) 
A Gagliardo-Nirenberg inequality yields 

\\A m h im+1) < H^^ll^^+^ll^lli 2 — 1>/2( ^ +1) +^-^^0 11^1^, (A.ll) 

m 

6 I am grateful to G. Ponce for pointing this result out to me. Note that the m = oo case is forbidden because an 
extra logi/3-term is needed [38 39J. It is this forbidden limit that ultimately prevents the closure of the gaps in the 
figures in Sf3] which allows trajectories to escape. 



3rd/09/09 (nsreg.tex) 



12 



which means that 

3/2m 



/ f \ 6 l Zm m + l I 

J m+ i ^ C ™\U |V(^ m )| 2 dVj J% m -^ 2m + L-*l m J^r \ ■ (A.12) 
With j3 m = m(m + 1), (|A.12|) can be used to form ft m +i 

L- 1 J | V(u m )\ 2 dV + L- 3 J m + vo 2m \ 



m+l 



( j — 3 t \l/2m ' 



\ 3/4/3 m 

(2m-l)/2(m+l) 



(A.13) 



which converts to 

c m (V 1 jf |V(u;-)| 2 + L- 3 J rn + Wq™^ > (^Y^ 3 • (A.14) 
This motivates us to re-write (IA.7P as 

-(L~ 3 j m ) < -^(l- 1 f \V(cj m )\ 2 dV + L- 3 J m + m 2 r 



+ C 2 , m (L Jm+l) m+1 ( Ji 



(m+l) 



+ c 3 , m tz7 L 3 J m + c 4 , 3 Wo m+1 + c 5, m ro o^m m X Gr . (A.15) 
Converting the J m into fi m and using Gr > 1 



< Mm { — f ^o +1 > ) + c 5,m (-77^-) fi I(m+l) + C 6 , m W Gr }(A.16) 

Using a Holder inequality on the central term on the right hand side (jA.16|) finally becomes 

f „ /Q \% 3_ 2(m+l) 1 

tim<n m { M-^l + c 2,«^o 2m "^i 2 ^n +c 3 ^oGr . (A.17) 

Cl, m V / 2( m+1 ) 



With no forcing the final term in (lA.17j) is proportional to w\. Converting to the dimensionless 





quantity D m = (w 1 Q m ) Qm already defined in (j2.ll) with a m = 2m/(4m — 3), finally gives 



{ \ (B ,\ 2m ( 4m + 1 )/ 3 \ 
{w Q a m )- 1 b m <D m { -^nM D 2 m + c 2:m D 2 n + c 3:m Gr} (A.18) 

[ Ci, m \ i^m / J 

with n = |(m + 1). ■ 
B Proof of Theorem Q] 

There exists a result of Foias, Guillope and Temam [32J, which uses higher derivatives. Define H n 

for n > 1 

«■„=/" |V n u| 2 dF, (B.l) 
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together with an integration of Leray's energy inequality 



T u ^ E 

m^L- 3 (H X ) T = {D X ) T < Gr 2 + -— -° . (B.2) 



T 

Then the result of Foias, Guillope and Temam [32J for n > 3 is 



H?- 1 ) < CnU^L- 1 

T 



Gr 2 + ^E 



T 

where Eq = E(Q) = i?o(0) is the initial energy. In the unforced case 



(B.3) 



H^) ^ . (B.4) 

A Sobolev inequality gives 

\\u\\ 2m <c\\V 2 u\\ a 2 \\u\\\- a (B.5) 

where a = 3(m — 1) /Am for m > 1. Moreover, the constant c can be taken as finite for each finite 
m because the m = oo case it is a bounded. Thus, taking n = 3 in (jB.3j) , which fixes the constant 
c rt , we have 

2m \ I / , \ m+3 \ 

| w ||^=s^ < c ^(# 1/5 ) 4(4m " 3) F* (4m - 3) ^ 

15(m-l) m+3 

< c (( J ff 3 1/5 ))^ (4m " 3) ( J H'i) 4(4m - 3) . (B.6) 

Using (fB~2l) and (fB~4|) this gives 

|u;|| 2 4 -- 3 ^ < [L- x Gr 2 + — £ J , (B.7) 

and thus the final result with an m-independent constant. In the unforced case 

There is also a way of reproducing the Gr 2 -estimate from Proposition [1] but with worse constants. 
Based on < JlJJJfii for n = \{m + 1), the relation in terms of the D n and D m is 

4m-3 1 



D l n < D^-'Dl™- 1 . (B.9) 

Inequality (1A.18I) is now divided by D 2 m & where 5 > ^m-i) ■ Noting that D m > 1 the Z^-term is 
handled as follows 

(2m-l)i5-l 



DlDt 2 ) T < ((^)^^(^V- - 



T 



(2m-l)<5-l\ / ,\ , 1 / 5 
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It follows that 



D m +1 \ 
D m J 



2m(4m+l)/3 



D 



< c Aym {D m ) + c b}rn {D\)+c^ m Gr + 0{T- v ) (B.ll 



where the coefficients from the Holder inequality have been absorbed into the constants. Define 

A m = 2m(4m + l)/3, and consider 



D 



m+l 



D 



m+l 



B r . 



D 



S/Ar, 



< 



m+l 



D°) + 



A \ V D 7 m / V A 

'- 1 m \ \ / / \ ^m 

where a Holder inequality has been used at the last step. The end result is 

D 5 m+1 ) < c 7 , m (D m ) + c 8 , m (D() + c 9 , m Gr + 0(T _1 ) . 



(B.12) 



(B.13) 



Because n = |(m + 1), when m = 1 then n = 1. Moreover, only when 6 = 1 does an estimate 
exist for {D\) through (IB. 21) . then (IB. 1311 is a generating inequality gives the Gr 2 -estimate but with 
worse constants. ■ 



C Proof of Theorem [2] 

With < fj, m < 1 we write 

/ nl— Mm 



< 



which becomes 



D„ 



D m +i 



1 — Mn 



D r 



Dm+l 



1— /in 



7-)l — Cm 

-^m+l 



m+l 



1 — yxn \ Mm 

Mm \ 



(D m +i) 



l-Mm 



> 



£)1— Mm \ \ 



r 



(D m+ \) T 



\ m / r j 



(C.l) 



(C.2) 



The estimate for the time average of (D m+ i) T from (|2.3p and the lower bound D m > 1 are now 
used to give 



D„ 



_D m +l 

This ends the proof of Theorem [2J 



\yav ' m 



1 — Mm \ 
Mm \ > . 

/ T 



(C.3) 
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